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Abstract 

We present results of a Monte Carlo simulation of a three dimensional Gross-Neveu 
model with SU(2) (g) SU(2) chiral symmetry at non-zero baryon chemical potential fi, 
corresponding to non-zero baryon density. For /i sufficiently large there is a sharp 
transition between a phase where the chiral symmetry is broken by a condensate (qq) 
and one where a scalar diquark condensate \{qq)\ 7^ 0. Global U(l) baryon number 
symmetry may remain unbroken, however, due to the absence of long range order 
in the phase of (qq). There is also tentative evidence for the formation of a weaker 
pseudoscalar diquark condensate in the high density phase, which violates parity. 
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1 Introduction 



The properties and behaviour of baryonic matter at high density have recently en- 
joyed renewed interest, with the observation that the ground state even at densities 
sufficiently high for chiral symmetry to be restored may be far from trivial |T| [0] . In 
brief, it is speculated that there is a regime of temperature and density for which the 
quarks have a large Fermi surface; kinetic considerations then suggest that if there 
is any attraction at all between quarks, then the free-particle vacuum is unstable 
with respect to condensation of diquark pairs from antipodal points on the surface, 
creating a gap between the Fermi energy and that of the lowest excited states. This 
is precisely a relativistic incarnation of the BCS instability in superconductors, with 
the diquarks playing the role of Cooper pairs. 

The condensation phenomenon can be modelled by assuming the interaction be- 
tween quarks is due to one-gluon exchange , or using effective four-fermion vertices 
resulting from the presence of instantons in the QCD vacuum 00. A recent cal- 
culation using a phenomenologically inspired Nambu - Jona-Lasinio (NJL) model is 
given in [Q. Estimates of the gap energy are of order lOOMeV, comparable with the 
chiral condensate at zero density and temperature. It is intrinsically interesting in 
a field theoretic sense because the diquark condensate may not be invariant under 
global U(l) rotations associated with baryon number, and hence its formation be an 
example of spontaneous breaking of a vectorlike symmetry. The condensate is also 
a gauge non-singlet, hence promoting a dynamical Higgs effect, and making some 
(or perhaps all 0) of the gluons massive, leading to "color superconductivity". The 
pattern of symmetry breaking is predicted to be critically sensitive to the number 
of light quark flavors present . Finally there may be interesting dynamical effects 
associated with the competition between chiral {qq) and diquark (gg) condensates as 
the current quark mass m and baryon chemical potential /i are varied. Phenomeno- 
logically, diquark condensation and its generalisations could have implications for the 
distribution of strangeness observed in heavy ion collisions 0, and perhaps in the 
suppression of neutrino emission from neutron star cores, resulting in slower cooling 
rates ii. 

So far calculations have relied on model approximations to the QCD interaction, 
with ad hoc form factors introduced, for instance, to model the effect of asymptotic 
freedom, which should suppress the condensate at large Fermi momentum, and hence 
large /i. Many assumptions about the nature and origin of the qq interaction are 
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necessary. This is a natural area, therefore, for numerical simulations of lattice QCD 
to make an impact. Recent simulations of QCD in a fixed gauge have found evidence 
for weak diquark binding at zero density [^. Unfortunately, with /i 7^ 0, the Euchdean 
QCD path integral measure detM(/i), where M is the fermion kinetic operator, is 
complex, making standard Monte Carlo simulation impossible. The most promising 
approach is to locate the zeros of the QCD grand canonical partition function in the 
complex fugacity plane avoiding the the problem of the complex action by generating 
the statistical ensemble using detM(jU = 0)e"'^'^ (where 5*^ is the gauge action) and 
reweighting with the factor ^ct M(At=o) achieve an overlap with the correct /i 7^ 
ensemble. This method, though exact in principle, appears extremely slow to converge 
in practice, and the results, though promising, are far from definitive P]. To date 
the lattice approach has only successfully been applied at non-zero density to toy 
field theories [3|T^[0. In this paper we present the first lattice study of diquark 
condensation in one such model, the Gross-Neveu (GN) model with SU(2) ® SU(2) 
chiral symmetry formulated in c? = 2 + 1 spacetime dimensions. 

Which features of the GN model make it worth studying? In brief |]I2| ||13[ : 



• For sufficiently strong coupling the model exhibits dynamical chiral symmetry 
breaking at zero temperature and density. 

• The spectrum of excitations contains both baryons (the elementary fermions), 
and mesons (composite fermion - anti-fermion states), which include a Gold- 
stone mode. 

• For 2 < (i < 4 the model has an interacting continuum limit. 

• When formulated on a lattice, the model has a real Euclidean action even for 
chemical potential 7^ 0, and hence can be simulated by standard Monte Carlo 
techniques. 

Therefore the model displays much of the essential physics except for color confine- 
ment. Even such a simple theory may be useful in addressing some of the issues 
raised above, including the behaviour of competing condensates, and the full phase 
diagram in the (/i, T) plane once temperature T > 0. Simulations at T = |]10| |jll] 
have revealed a first order chiral symmetry restoring transition at /ic — ^f-, where 
is the physical (ie. constituent) fermion mass at = 0; this is in marked contrast 
to the pathological behaviour observed in approximations to QCD when the measure 



is held real, such as the quenched approximation, where (ic ~ "^7r/2 . 
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In the next section we will describe the lattice model in detail, and identify which 
diquark operators we have studied using both lattice and continuum notation. Our 
simulations and results are described in section 3, and a short discussion given in 
section 4. 



2 The Lattice Model 
2.1 Action and Symmetries 

The lattice model we have simulated has the following Euclidean path integral: 

Z = J DaDir det(M^M[a, vf]) exp (^-^(^^ + ^•^) j ' (1) 

where a and the triplet tt are real auxiliary fields defined on the dual sites x of a three 
dimensional Euclidean lattice. The determinant factor can be viewed as having been 
obtained by integrating over Grassmann variables x? C? C with the following action: 



with 
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Sf,r = xMW,n]x + CM*[a,n]C 



u=l,2 



+ k,EK^)^"'' + ^^(^)^(5^)-^"'^)- (3) 

O ,- , 

{x,x) 

The parameters are the bare fermion mass m, the chemical potential /i and the 
coupling g^. The r are Pauli matrices with indices p,q = 1,2. The symbols ?7jy(x) 
denote the phases (— \-x^-i^ ^^^^ ^j-^g phase (^—i) ^'0+^1+^2^ g^^^j the set of 8 

dual lattice sites neighbouring x. In this form we see that x ^"^^ C can be identified 
with isodoublet staggered lattice fermion fields defined on the sites x. The integration 
measure in (|I|) is manifestly real|^; note that x and ( have opposite-signed couplings 
to the fields tti and 713. It is straightforward to integrate over the auxiliary fields 
a and vr to recover an action which contains fermions which self-interact via a four 
point interaction oc g^; more details are given for the corresponding four- dimensional 
model in 



^Note also that detM = detT2M*T2 is real m§ 
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In the limit m — > 0, the model is invariant under a global symmetry akin to the 
SU(2)^ ® SU(2)^ of the continuum NJL model. Defining projection operators onto 
even and odd sublattices Ve/o by 

^e/o(x) = i(l±£(x)), (4) 



we have 



C^{VeV + VoU)C ; C^aVeU^ + VoV^) (5) 
$ = (al + ivr.r) ^ V^U'' 

with U,V E SU(2). Now for T = yU = 0, a mean field treatment (which for a model 
with fiavors of x ^iid C, is equivalent to the leading order of an expansion in 1 / A^) 
shows that this SU(2) ® SU(2) symmetry is spontaneously broken to 811(2)^^^^^^^ 
by the generation of a condensate Sq = (cr) = p-(xx) and physical fermion mass 
TTLf = So, the symmetry being broken for coupling > (7^ ~ 1.0. Hence Sq defines a 
physical scale in cutoff units: a continuum limit exists as from either phase. 

Quantum corrections to this picture can be calculated as higher order terms in the 
expansion, which is renormalisable about g^ = g1 [|l^[|l3l. The numerical results 



in this paper have been obtained for = 1 and a value g = 2.0, corresponding to a 
zero density theory deep in the broken phase, with Sq = 0.706(1) in units of inverse 
lattice spacing, ie. rather far from the continuum limit. 

We can also identify two other symmetries of (|1|J5), both of which hold even for 
m 7^ 0. First, there is the global U(l) of baryon number: 

X,C^e*"x,C ; xX^e-'^xX- (6) 

The chemical potential ^ couples to the conserved charge associated with this sym- 
metry. Next, there is a discrete parity symmetry appropriate to 2+1 dimensions: 

a; = (a;o,a;i,a;2) ^ x' = (xq, 1 - xi, 0:2) 

X, C(x) ^ (-i)<+^^x, a^') ; X, Cix) ^ {-ir'^^-'-x, Cix') (7) 

(t{x) h- * cr(x') ; 7f(x) 1-^ — 7r(a;'), 

whence the identification of a as scalar and tt as pseudoscalar. 

For T, /X 7^ a mean field solution is also known |T^|T^, in which S(/i,T) is 



expressed in terms of Eq. At zero temperature the basic feature is that E remains 
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constant as /i is increased up to a critical value /Xc = ^O; whereupon S falls sharply to 
zero (in the chiral limit), signalling a first order chiral symmetry restoring transition. 
Monte Carlo simulations [|I0] support this picture, even when massless Goldstone 
excitations (with the quantum numbers of the vr field) are present, and N takes the 



minimal value = 1 11 



The reality of the measure in (|I|) is an artifact of our introducing conjugate flavors 
X and (. It is worth discussing how the Monte Carlo simulation manages to reproduce 
the main features of the expected behaviour of the model once 7^ 0, while similar 
simulations of QCD with quarks and conjugate quarks fail dramatically |T^[|r^. In 



brief, it is because interactions between x and ( are negligible so that a light bound x( 
state, which would carry baryon number, does not form. The strongest interactions 
are in the singlet channels XX and ((, which are dominated by disconnected "bubble" 
diagrams in the 1/A^ expansion, and it is in these channels that Goldstone poles form 
||19| . Similarly, one might wonder how a diquark condensate, which is not invariant 
under the vectorlike U(l) global symmetry of baryon number, could ever form in a 
model with a real measure, since such a breaking is usually forbidden by the Vafa- 
Witten theorem |2^. The resolution is that the theorem does not hold for models 
with Yukawa couplings to scalar degrees of freedom, such as the coupling to a in (|^). 

2.2 Constructing Diquarks 

Next we discuss the operators used to form diquark pairs. We begin by expressing 
the operators in terms of the lattice fields x (we did not perform measurements in the 
( sector, but the results here are identical). In this study we have used four different 
local diquark operators: 



NSS 



: xnx)xnx) NSP : £(x)x^(x)xn^) 
''(x)rf x'(x) SP : e{x)x''{x)Trx'{x). 



SS: X 

Operators NSS and NSP are symmetric in both spatial and isospin indices, and so 
would only form a non-zero condensate if extra flavors are introduced, ie. > 1. In 
this case the operator could be written, with the indices i,j running from 1 to N, 
^Xi^ijXj- For N = 1 this vanishes identically, but as described below, the diagram 
which would contribute for N > 2 can also be measured on an ensemble generated 
with = 1 in a generalisation of the quenched approximation; we will refer to such 
operators as non-spectral diquarks. Operators SS and SP, on the other hand, are 
symmetric in space but antisymmetric in isospin, and hence can form even for N = 1; 
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we will refer to these as spectral diquarks. States NSS and SS are even under parity 
(0), and so will be called scalar, while NSP and SP are odd, and hence pseudoscalar. 
Note that NSS and NSP are not invariant under the chiral rotations or indeed 
under the remnant isospin symmetry in the chirally broken phase, and all four diquark 
operators not invariant under the U(l) of baryon number (j^). 

To gain further insight it is useful to transform to a basis in which continuum-like 
spinor indices are shown. The recipe for staggered fermions in three dimensions is 



well-known First we make a unitary transformation to fields u and d defined on 



a lattice of twice the spacing of the original, with site labels y 

1 ^„ . 1 



(?/) = Z7^ErTx(^;i/) ; ^"%) = ^E^rx(^;?/)- (9) 



The indices a, a each run from 1 to 2. Here A is a lattice vector with entries either 
or 1, so that each site x on the original lattice corresponds to a unique combination 
of y and A. The 2x2 matrices F and B are defined by 

TA = r^'r,^'Ti- ■ Ba = {-t,)^\-t,)^K-t,)^'. (10) 

Now, with the definition 

we see that q may be interpreted as a four-component spinor operated on by Dirac 
matrices defined by: 

with two flavors counted by the latin index a; this extra flavor degree of freedom is 
due to the species doubling inherent in the lattice approach. It is now possible to 
recast the fermion matrix (^ including the Yukawa interactions in continuum- like 
notation [^[|T^]. Here we will content ouselves with writing down expressions for the 



diquark operators (P) in the g-basis: 

NSS : -ig^(C75 ® ra ® l)q NSP : -ig^(C ® Ts ® l)q 
SS : -ig^(C75 T2)q SP : ~iq^ {C T2 Tg)?. (13) 

Here the first matrix in the direct product acts on the 4-spinor indices, the second 
on the implicit flavor space indexed by a, and the third on the explicit isospin degree 
of freedom indexed by p, q. The matrix C is the anti-unitary charge conjugation 
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matrix defined by C'j^C'^ = —7^; in our explicit basis (|T2D C = 7072- Denoting a 
symmetric state by s and an antisymmetric by a, we see operators NSS, NSP are 
a (S> a ® s = s, and SS and SP a ® a ® a = a, corroborating our identification of 
spectral and non-spectral operators made earlier. 

We can also rewrite the symmetries in the g-basis. The SU(2) ® SU(2) chiral 
symmetry (j^) becomes 

QL^UqL ; Qr^ Vqr with Ql/r = ^(1 ± (75 (g) 1 (g) 1)) q, (14) 

the U(l) of baryon number (|^) 

g exp(za(l (g> 1 (g) l))g ; g gexp(-m(l (g) 1 (g 1)), (15) 

while the parity transformation (|^) is now 

g(a;) (7175 (g 1 ® l)g(x') ; g(a;) g(a;')(757i ® 1 ® !)• (16) 

The parity assignments of the diquark operators (|T3|) can be checked by observing 
that 7i75 is a symmetric matrix. 

3 Simulations 

We have performed simulations of the model (|l],0) on x Lt lattices with bare mass 
m = 0.01 at a coupling = 0.5, with Ls varying between 8 and 24, and Lt between 
24 and 40. The chemical potential /i was varied between and 1.2, with a critical 
value for chiral symmetry restoration ~ 0.65. The simulation method is a standard 
hybrid Monte Carlo algorithm. We monitored the chiral condensate (xx)) the 
baryon density n defined by 

n = ^ (x^(x)e'^x^(a; + 0) + x^(x)e~V(a; - 0)) ~ (g(7o ® 1 ® l)q) (17) 

using stochastic estimators. 

To estimate diquark condensates we chose an indirect method: the two point 
function (gg(O)gg(x)) was measured as the expectation of the product of two fermion 
propagators as in QCD baryon spectroscopy, and the condensate extracted by assum- 
ing clustering at large spacetime separation, ie. 

{qq{0)qq{x)) = {qq{0)qq{x))c + {qq){qq), (18) 
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Figure 1: Spectral scalar correlators for two values of /j, in the broken phase on 16^ x 24 
lattices. 

the latter term on the right being proportional to \{qq)\'^- A non-zero condensate 
signal therefore shows up as a plateau in the timeslice correlator G{t) at large time 
separation, and may be extracted by a fit of the form 

Git) = J2 {qqmq{t,Xi,X2)) = Aexp{-M+t)+Bexp{-M_{Lt-t))+C\L,)\{qq)\', 

Xl,X2 

(19) 

with M± the masses of forward and backward propagating diquark states respectively. 
The value of the constant C{Ls) is not determined a priori. Naively one expects 
C = Ls, and hence the plateau height to be an extensive quantity. 

We found by measuring the chiral condensate that chiral symmetry is restored for 
chemical potential greater than fi^ — 0.65. The lattice virtually saturates with two x 
quarks per site at fi = 1.5 corresponding to a number density n = 2. 

The spectral scalar propagators in the broken phase n < are shown in Fig. |I[ 
For /i = 0.0 the spectral scalar condensate signal is vanishing while for /i = 0.6 it is 
very small but non- vanishing. We found no signal for the spectral pseudoscalar in the 
broken phase. The symmetric phase propagators and their fits are shown in Fig. |[ 
Noting the change in scale on the G{t) axis and comparing with Fig. |I| we see that 
in the chirally symmetric phase a strong signal develops for the spectral scalar which 
increases with chemical potential. The spectral pseudoscalar signal is weaker but 
distinctly non-zero and again increases with /x. The non-spectral propagators are not 
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Figure 2: Spectral scalar and pseudoscalar correlators and their fits for two /x values 
in the chirally symmetric phase on 16^ x 24 lattices. 



shown here although both the scalar and pseudoscalar closely follow the behaviour of 
the spectral pseudoscalar. 

The bulk of our results are from 16^ x 24 lattices. However, to assess possible 
systematic effects arising from fitting the diquark propagators on lattices of finite 
temporal lattice extent Lt we repeated simulations at = 0.4 on 16^ x 32, and 
/i = 0.8 on 16^ X 32 and 16^ x 40 lattices. The results at = 0.8 are shown in Fig. 
^, together with results obtained for free fermions (ie. with g"^ = 0) on the same 
lattices with the same chemical potential. The time axis has been normalised by 
L^^ to aid comparison. Note that the plateau height in the middle of the lattice is 
almost constant as Lt increases, consistent with | (qq) \ ^ 0, whereas in the free case 
the signal falls monotonically. We found that the fits to the spectral scalar correlator 
in the symmetric phase (/i = 0.8) were very stable, yielding estimates for C|(gg|) 
of 0.425(1) for U = 24, 0.424(1) for Lt = 32, and 0.436(1) for Lt = 40. This is 
strong evidence for a non-vanishing condensate. The smaller parity violating signal, 
on the other hand, decreased by about 25% from 0.197(1) at Lt = 24, via 0.167(1) 
at Lt = 32, to 0.155(1) at Lt = 40. In the broken phase (/i = 0.4) the spectral scalar 
C\{qq)\ was 0.065(1) for Lt = 24 but consistent with zero on the Lt = 32 lattice 
suggesting that the \{qq)\ ^ result in the broken phase is at least in part due to 
finite lattice size. 

The trend in the diquark condensates as we pass from broken to symmetric phase 
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Figure 3: Spectral scalar correlators measured on 16^ x Lt lattices at /i = 0.8, together 
with the corresponding results for free fermions. Note that the t axis has been rescaled. 

is clear from Fig. ^. There is a critical value, Hc — 0.65, at which the chiral (qq) 
condensate falls sharply and the number density begins to rise from zero. The spectral 
scalar C\{qq)\ condensate rises slowly from zero for 0.4 < fi < fic, jumps upwards 
discontinuously close to fic, and continues to rise steadily. For /i > 1.0 as the number 
density approached saturation the scalar qq propagators became distorted and we 
were unable to achieve satisfactory fits. 

The existence of a non-zero (qq) condensate in the broken phase could in principle 
convey physical information about the onset of a nucleon liquid type phase. A non- 
zero condensate implies that there is a Fermi surface and therefore non-zero number 
density. However, from the study with varying Lt discussed above, it is clear that 
further study is required to determine whether there is any non-zero signal associated 
with the onset of a nucleon liquid phase. 

The spectral pseudoscalar C\{qq)\ signal was consistent with zero for fi < fi^ 
and no fits to the propagators were possible but for n > this condensate is non- 
zero and increases with /i. It is considerably smaller in magnitude than the spectral 
scalar. This pseudoscalar condensate is parity violating, and therefore it would be 
remarkable if the signal in the chirally symmetric phase were to remain non-zero in 
the large volume limit. 

To determine the behaviour with spatial volume we next performed a series of 
runs on x 24 lattices at /i = 0.8, with Ls varying from 8 to 24. To our surprise 
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Figure 4: Overview of the observables (xx); ^^^^ SS and SP diquark condensate 
signals C\{qq) \ as functions of fi, obtained on a 16^ x 24 lattice, showing the onset of 
diquark condensation at fic — 0.65. 




Figure 5: Spectral scalar diquark correlator at /i = 0.8 for both free and interacting 
quarks on x 24 lattices. 



12 



0.80 



0.60 



S 0.40 



0.20 



0.00 



I spectral scalar qq 
' M spectral scalar qq 



0.30 



0.90 



0.50 0.70 
chemical potential 

Figure 6: Masses of forward and backward moving states and M_ obtained from 
fits to the spectral scalar propagators on a 16^ x 24 lattice plotted as a function of fi. 



we found very little change as Ls increased, the fitted value of C\{qq)\ more or less 
saturating for Lg > 16. Results for Lg = 8, 16, 24 are shown in Fig. ^, together with 
the fit for Lg = 24. For comparison we also plot the equivalent correlators for free 
fermions at yU = 0.8; it is striking that for this case the trend as Lg increases is in the 
opposite direction. We conclude that the constant C is roughly independent of L^; 
this puzzling behaviour will be further discussed in the final section. 

The diquark masses M+ and M_ obtained from the fits to the spectral scalar 
qq propagator are plotted as a function of /i in Fig. ^. We found that the forward 
propagating M+ states were more difficult to extract from the fits than the backward 
propagating M_ states, as refiected in the size of the error bars. The trend in the 
data is clear with small masses in the symmetric phase and large masses in the 
broken phase. The approximate linear decrease of M_ with fi refiects a similar trend 
observed in the physical fermion mass in previous simulations |TT|. Notice that 
M_ is approximately constant for /i > fic- The observed qq states have masses around 
0.3 which is light in comparison to the zero density fermion mass rrif = Sq — 0.7, 
which defines the ratio of physical to cutoff scales. 

The spectral pseudoscalar and both nonspectral C\{qq)\ condensates are plotted 
as a function of fi in Fig. |^. The signal for these three condensates is much weaker 
than the signal for the spectral scalar condensate. The parity-violating spectral pseu- 
doscalar and the two non-spectral condensates are almost identical in magnitude and 
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Figure 7: The three weaker diquark condensates C\{qq) \ plotted in the high density 
phase as a function of fi. 

have remarkably similar behaviour as a function of /i. We also remark that we have 
observed no evidence for spontaneous violation of isospin on inspection of the fermion 
propagators. 

To determine the effect of the bare quark mass on the (qq) condensates we per- 
formed simulations on 16^ x 24 lattices for m = 0.02,0.05,0.08,0.1 at = 0.8. The 
values of the condensates were consistent within errors with the m = 0.01 result 
C|(gg) 1=0.424(7) for m = 0.02,0.05 but the signal decreased by around 20% to 
0.342(10) for m = 0.08 and by a further 20% to 0.272(6) for m = 0.1. In fact 
all of the diquark condensates decreased as we moved further from the chiral limit. 

4 Summary and Outlook 

We have simulated the 3d GN model with SU(2)®SU(2) chiral symmetry at non-zero 
chemical potential and found evidence suggestive of a diquark condensation, namely 
that the two-point correlation function exhibits clustering at large temporal separa- 
tion. In order to quantify the measurement, however, it will be necessary to have a 
numerical estimate for the constant C in (|19|) . We have up to now ignored the pos- 
sible manifestations of fluctuations in the phase of the condensate. The spontaneous 
breaking of the U(l) symmetry (|]) usually implies the existence of a Goldstone mode 
associated with long wavelength fluctuations in the phase of the (qq) condensate; this 



14 



is not accessible by the methods presented here, since G{t) is real. However, in the ab- 
sence of an explicit diquark source (analogous to a bare mass for the case of (qq)), one 
would expect large finite volume effects, generically proportional to L~^'^'~'^\ where d' 
is the dimension of the effective field theory describing the fluctuations of the order 
parameter We therefore speculate that the effective field theory describing the 
spatial correlations of the diquark correlator has d' = d — 1 = 2, and that there is no 



long range order in the spatial directions . A similar distinction between temporal 
and spatial correlations has been observed in instanton liquid models . This would 
account for the volume-independence of C. Strictly speaking, therefore, we have not 
observed a true condensation. This interpretation of our results can be tested by 
simulations of the equivalent 3 + 1 dimensional model [|l^], where we would expect 
long range order. 

In a BCS condensation, only quark pairs with momenta ±p close to the Fermi 
surface contribute to the condensate. Therefore we expect that C^, which must 
depend on the number of participating qq states, will be proportional to the area of 
the Fermi surface, ie. oc fi'^~'^. Note that the spectral scalar data in our work is 
well fitted by C\{qq) \ oc y/fi for 0.75 < < 9; fits over a wider range of fi, obtained 
from simulations closer to the continuum limit, would help to confirm this behaviour. 
If we assume that C ~ 0{y/fi), we can claim to have detected a strong signal for a 
spectral scalar | (qq) \ condensate in the chirally symmetric phase which is of the same 
order of magnitude as the chiral (qq) condensate in the broken phase. This would 
imply that the energy gap associated with diquark condensation is 0(So) and hence 
support the predictions of . 

We have observed a signal for | (qq) | 7^ in the broken phase possibly associated 
with the onset of a nucleon liquid phase expected to occur for /i ^ fi^, but further 
study of finite size effects is needed for confirmation. We have also detected a weak 
signal for a parity violating | (qq) \ condensate in the chirally symmetric phase. This 
may be a three dimensional analogue of the axial diquark condensate reported in [Q. 
The remarks about finite size systematics apply here too; however this fascinating 
phenomenon could also be investigated by checking for the presence or absence of 
parity doubling in the spectrum of spin-1 states. In a parity-symmetric vacuum we 
expect odd and even parity spin-1 states to be degenerate0. This can be observed in 
the spectroscopy of the 3d Thirring model [p^ . 

The two-point method used here to measure the {qq) condensates requires the 
^ This was pointed out to us by M. Teper. 
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diquark operators (^) to be local in the x fields; this excludes potentially interesting 
operators in this and related models. The introduction of an explicit diquark source to 
study the one-point function directly is thus desirable, both to calibrate our indirect 
two point function measurements, and to extend our study to other diquark operators 
and the T 7^ regime. 
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